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La presente invention conceme une methode pour former plus rapidement un 
modele numerique stochastique de type Gaussien ou apparente, representatif d'un milieu 
heterogene poreux (tel qu*un gisement d'hydrocarbures par exemple) cale par rapport a des 
donnees dites dynamiques caracteristiques du deplacement des fluides dans )e milieu telles 
que par exemple des donnees de production (pressions obtenues a partir d'essais de puits, 
debits, etc.). 

La methode selon I'invention trouve des applications dans la modelisation de zones 

souterraines ou il s'agit de generer des representations montrant comment est distribuee 

■ I ' 

une certaine grandeur physique dans une zone du sous-sol (la permeabilite notafhment), 

compatibles au mieux avec des donnees observees ou mesurees, dans le but par exemple 

d*en favoriser Texploitation. • - 

ETAT DE LA TECHNIQUE 

L' optimisation dans un contexte stochastique consiste a deteiTnirier des realisations 
d'un modele stochastique qui satisfont un ensemble de donnees observees sur le terrain. En 
ingenierie de reservoir, les realisations a identifier correspondent a des representations, 
dans le champ reservoir, de la distribution de proprietes de transport telles que la 
permeability ou porosite. Ces realisations forment des modeles numeriques de reserx'oir. 
Les donnees disponibles sont, par exemple, des mesures ponctuelles de permeabilite ou 
porosite, un modele de variabilite spatiale determine selon des mesures ponctuelles ou 
encore des donnees directement liees aux ecoulements des fluides dans un r6sei*voir 
soulerrain, c'est a dire des pressions, des temps de percee, des debits, etc. Ces demieres 
sont souvent non lineairement reliees aux proprietes physiques a modejiser. Une realisation 
tiree au hasard n'est pas en general en adequation avec Tensemble des donnees collectees. 
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La coherence vis a vis des donnees est integree dans le modele par le biais d'une procedure 
inverse : 

- Tarantola, A., « Inverse problem theory - Methods for data fitting and model parameter 
estimation », Elsevier Science Publishers, 1987, 

5 Pour ce faire, la technique la plus simple est celle de Tessai et de Terreur. Cette 

approche consiste a tirer au hasard des realisations jusqu'a I'obtention d'une realisation 
satisfaisant aux donnees. Elle offre I'avantage de la conservation du modele de variabilite 
spatiale, mais requiert un temps calcul prohibitif. Aussi n'est-elle que tres rarement utilisee 
en pratique. 

10 Une alternative souvent preferee prend appui sur le calcul des gradients. Les 

methodes de gradients pemiettent de modifier une realisation initiale suivant une direction 
de recherche, cette demiere etant estim6e a partir des gradients. Les modifications sont 

appoitees iterativement jusqu'a ce que le calage des donnees soit juge acceptable. Les 
methodes de gradients sont attrayantes de par leur efficacite. Neanmoins, elles deviennent 
15 inadaptees des lors que le nombre de parametres, c'est-a-dire le nombre de valeurs de 
permeabilite et de porosite constituant le modele numerique, est important. En outre, elles 
ne permettent pas de modifier les realisations en respectant la structure spatiale du modele 
stochastique. 

Plus recemment, une technique de parametrage geostatistique, a ete introduite pour 
20 contraindre, pai^ deformation graduelle, les realisations stochastiques a des donnees dont 
elles dependent de maniere non lineaire. Elle a fait I'objet des brevets FR 2.780.798 et 
FR2.795.841 du demandeur, et des publications suivantes, notamment : 

- Hu, L.Y., 2000, Gradual deformation and iterative calibration of Gaussian-related 
stochastic models: Math, Geology, VoL32, No.l. 

25 - Le Ravalec, M. et al. 2000, The EFT moving average (FFT-MA) generator: An 
efficient numerical method for generating and conditioning Gaussian simulations: 
Math. Geology, VoL32, No.6. 

- * Hu, L.Y., Blanc, G. and Noetinger, B. (2001): Gradual deformation and iterative 

calibration of sequential stochastic simulations. Math. Geology, Vol. 33, No.4. 
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Cette methode a ete appliquee avec succes a divers cas notamment a partir de 
donnees provenant de champs d'exploitation petroliere, comme decrit dans les documents 
suivants : 

Roggero, F, et al,1998, Gradual deformation of continuous geostatistical models for 
5 history matching, Paper SPE 49004: Proc. SPE Annual Technical Conference and 

Exhibition, New Orleans. 

Hu, L.Y. et al. 1998, Constraining a reservoir facies model to dynamic data using a 
gradual deformation method, Paper B-01: Proc. 6th European Conference on 
Mathematics of Oil Recovery (ECMOR VI), 8-11 September 1998, Peebles, Scotland. 

10 Comme on va le rappeler plus en d6tail ci-apres, la methode de deformation 

graduelle permet de modifier graduellement line realisation d'un modele stochastique de 
type Gaussien ou de type appaiente au modele Gaussien tout en respectant la staicture 
spatiale du modele. 

Optimisation stochastique 
15 Soit A'^' Af) I'ensemble des donnees sur le teniin et 

■ ^ 

f - (fufi reponses correspondantes simulees pour une realisation, z d'un 

modele stochastique Z donne. En general, les reponses / = (/, /2 /m ) sont obtenues en 

resolvant numeriquement le probleme direct. Ainsi, si z represente un champ de 
permeabilite, les donnees /peu vent etre des mesures de pression. Dans ce cas, elles sont 
20 simulees a paitir d'un simulateur d'ecoulement. L'objectif d'une optimisation stochastique 
est de produire des realisations de Z qui reduisent les differences entre les donnees 
observdes et les reponses correspondantes simulees numeriquement. Ces differences sont 
mesurees par la fonction objectif suivante : 

25 Les coefficients (i>tn sont des poids attribues aux donnees f^^ . Les sont des 

fonctions de la realisation z. En ce sens, la minimisation de la fonction objectif est un 
probleme a plusieurs variables. 
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Soit N ]q nombre de mailles ou de composantes formant la realisation z. N est 
souvent tres grand (lO'^ --10'^). II est done tres difficile de mener une optimisation 
directement par rapport aux composantes de z. De plus, la realisation z, meme modifiee, 
doit rester une realisation de Z. La methode de deformation graduelle permet de 
5 s'affranchir de ces difficultes. 

Recherche aleatoire d partir de la methode de deformation gi'aduelle 

On considere a present un champ aleatoire Z qui peut etre transforme en champ 
aleatoire Gaussien Y. La technique de deformation graduelle permet de construire une 
chaine continue de realisations en combinant une realisation initiale 3^0 de Y avec une autre 
10 realisation wi, dite complementaire, de 7, ui etant independante de yo (figure la). Les 
coefficients de combinaison sont par exemple cos(0 et sin(0 et la realisation combinee 
verifie la relation : 

y{t)~ y^cost ^u^smt 

ou t est le parametre de defomiation. 

15 Une autre technique de construction de chaines de realisations consiste a combiner 

la realisation initiale avec non pas une, mais P realisations complementaires Up (p=l,P) de 
Y (figure lb), Les coefficients de la combinaison sont tels que la somme de leurs carres est 
egale a 1. 

Des lors que la chaine est elaboree, on peut rexplorer en variant le parametre de 
20 deformation t et tenter d'identifier parmi toutes les realisations de cette chaine celle qui 
minimise la fonction objectif. Cette minimisation se fait par rapport a /. Le parametrage 
suivant la methode de deformation graduelle permet de reduire le nombre de dimensions 
du problenie de N a 1, ou //est le nombre de valeurs constituant le champ a contraindre. 
De plus, la somme des coefficients de combinaison au caixe etant 1, la realisation 
25 optimisee y{topt) est encore une realisation de Y : elle suit le meme modele de variabilite 
spatiale que toutes les realisations de Y. 

Cependant, si Ton restreint Texploration de Tespace des realisations a une unique 
chaine, on limite severement nos possibilites de. diminuer suffisamment la fonction 
objectif, II faut done repeter la procedure decrite d-dessus, mais avec de nouvelles chaines 
30 de realisations. Ces chaines de realisations sont construites successivement en combinant 



1 er depot 

5 



une realisation initiale qui est ici la realisation optimale determinee a Titeration precedente, 
avec une realisation complementaire de a chaque fois tiree au hasard. Ainsi, a Titeration 
/, la chaine continue de realisation s*ecrit : 



independantes de Y. Ces dernieres sont aussi independantes de vq. Cette foiTnulation 
implique que la chaine de realisations correspond a une hyper ellipse de dimension N. 

En minimisant la fonction objectif par rapport a /, on amel.iore, ou au moins 
preserve, le calage des donnees a chaque fois qu'une nouvelle chaine de realisations est 

10 exploree. Cette procedure de recherche de minimum iterative est poursuivie tant que le 
calage des donnees n'est pas satisfaisant, Le cote aleatoire de la methode tient au fait qu'a 
chaque iteration, la realisation complementaire est tiree au hasard. De fait, la direction de 
recherche dans laquelle on part depuis la realisation optimisee a Tetape precedente est 
aleatoire. En effet, la direction de recherche, pour une chaine donnee et a partir de la 

15 realisation optimale definie precedemment, est : 



Cette direction de recherche depend uniquement de w/. De plus, U[ etant 
independante des realisations complementaires deja generees «, aussi de yo^ la 

direction de recherche au depart de chaque nouvelle chaine est orthogonale a la tangente 
20 definie pour la chaine precedente au meme point (figure 2). Bien qu'il puisse paraitre 
adequat de selectionner une direction de recherche orthogonale a cette tangente, il existe 
une infinite de directions orthogonales possibles. 

L'experience montre que, apres quelques iterations, les nouvelles directions de 
recherche ne contribuent plus significativement h la d6croissance de la fonction objectif 
25 (figure 6), 

II a ete envisage aussi de combiner la realisation initiale non pas a une, mais a 
plusieurs realisations complementaires. Dans ce cas, le nombre de parametres de 
defonnation augmente : il egale le nombre de realisations complementaires mises en jeu 



yi (0 = 3'/~i ^os/ + III ^ • 
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yi-i est la realisation optimale definie a Fiteration M et les ui sont des realisations 
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lors d'une combinaison gi-aduelle. Bien que le processus d' optimisation sort alors plus 
flexible, il faut pouvoir gerer plusieurs parametres ce qui n'est pas forcement evident. Par 
ailleurs, un tel processus n'est pas necessairement plus efficace, car il peut dependre de 
I'execution d'un plus grand nombre de simulations directes d'ecoulement. 

La methode selon rinvention 

La m6thode selon I'invention permet de former plus rapidement un modele 
numerique representatif de la distribution d'une grandeur physique dans un milieu 
heterogene poreux, tel qu'une zone souterraine (reservoirs petroliers, aquiferes, etc.), 
contraint par des donnees collectees dans le milieu (donnees dynamiques caracteristiques 
du deplacement des fluides dans le milieu, collectees par des mesures (dans des puits de 
production, d'injection ou d'observation par exemple) ou des observations prealables. 

Elle comporte un processus iteratif de deformation graduelle ou Ton combine 
lineairement h une realisation initiale d'une partie au moins du modele choisi du milieu, au 
moins une deuxi^me realisation independante de la realisation initiale, les coefficients de 
cette combinaison lineaire etant tels que la somme de leurs carres est egale a L et I'on 
minimise une fonction objectif mesurant I'ecart entre un jeu de donnees non lineaires 
deduites de la dite combinaison au moyen d'un simulateur du milieu, et les dites donnees 
geologiques et dynamiques par ajustement des coefficients de la combinaison, le processus 
iteratif etant repete jusqu'S obtenir une realisation optimale du modele stochastique. 

La methode est caracterisee en ce que Ton accelere la vitesse de deformation 
graduelle vers le modele optimal representatif du milieu, en selectionnant comme 
deuxieme realisation a combiner avec la realisation initiale, au moins une realisation 
composite obtenue en selectionnant au prealable une direction de descente definie en 
fonction des gradients de la fonction objectif par rapport a toutes les composantes de la dite 
realisation initiale. 

On obtient la realisation composite par exemple par combinaison lineaire d'un jeu 
de realisations independantes du modele, les coefficients de la combinaison etant calcules 
de fa^on que la direction de descente a paitir de la realisation initiale y soit le plus proche 
possible de celle definie par les gradients de la fonction objectif par rapport h toutes les 
composantes de la realisation initiale. 
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On conduit T optimisation par exemple a paitir d'un parametre de deformation, qui 
controle la combinaison entre la realisation initiale et la realisation composite. 

Dans le cas ou la dite combinaison n'affecte qu'une partie de la realisation initiale, 
on applique le processus iteratif de defonnation graduelle a un bruit blanc Gaussien utilise 
pour generer une realisation gaussienne et on determine les derivees de la fonction objectif 
par rapport aux composantes du bruit blanc Gaussien. 

Suivant un mode de mise en oeuvre, on combine la realisation initiale avec un 
certain nombre M de realisations composites, toutes obtenues par composition a partir de 
F„, realisations independantes de 7, 1'optimisation faisant intervenir M parametres. 

En d'autres termes, la methode comprend essentiellement un nouveau schema de 
combinaison graduelle qui tient compte de I'information apportee par les gradients au point 
initial de toute chame de realisations. Pour construire une chaine, on part toujours d'une 
realisation de depart et d'un ensemble de realisations complementaires, toutes 
independantes et issues du meme modele stochastique. Par contre, on ne combine pas 
directement la realisation de depart avec les realisations complementaires. Les realisations 
complementaires offrent la possibilite d'explorer Tespace des realisations dans des 
directions differentes. Ces directions ne sont pas equivalentes : certaines pemiettent de se 
dinger davantage vers F optimum. A ce stade, on elabore une realisation, dite composite, 
en combinant les realisations complementaires uniquement. Puis, on cree une chaine de 
realisations a partir de la realisation de depart (initiale) et de cette realisation composite. 
Cette chaine, comme celle qui etait propos€e dans le cas de base de la deformation 
graduelle, peut etre exploree a I'aide d'un seul parametre de deformation. 

La realisation composite est construite de fa5on a proposer une direction de 
recherche le long de la chaine aussi proche que possible de la direction de descente donnee 
par les gradients. Comme mentionne plus haut, toutes les realisations complementaires ne 
sont pas equivalentes : la realisation composite prend le meilleur de chaque realisation 
complementaire. 

La methode permet d'arriver plus rapidement a la formation du modele numerique 
representatif du milieu. 
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Presentation des figures 

D'autres caracteristiques et avantages de la methode selon I'invention apparaitront 
plus precisement a la lecture de la description ci-apres d'un exemple non limitatif 
d' application, en se referant aux dessins annexes ou : 

- les figures la, lb montrent des schemas de deformation graduelle (dits GDM) deja 
' conn us ; 

- la figurelc montre le schema de deformation graduelle (dit GBC) coirespondant h la 
methode selon I'invention ; 

- la figure 2 represente des chaines de realisations dans un espace euclidien a N 
dimensions, ou la tangente au niveau de la realisation optimisee pour une chajne de 
realisations l-l (/?C;., ) est orthogonale a la direction de recherche pour la realisation de 
depart de la chaine de realisations suivante / ( RC, ); 

la figure 3 montre la projection de la direction de descente v dans le sous-espace U 
defini par la base orthonormale formee par P realisations independantes (u^ Up ) ; 

les figures 4A k 4C montrent des exemples de distributions de permgabilite 
respectivement pour les realisations de reference, initiale et contrainte aux donn^es de 
pression; 

les figure 5A a 5E montrent respectivement les variations en fonction du temps des 
pressions simuldes respectivement aux cinq puits (BHP-OBSl, BHP-OBS2, BHP- 
PROl, BHP-OBS3, BHP-OBS4) des figures 4, respectivement pour les distributions de 
permeabilite de ref6rence (data), initiale et contrainte (match) ; et 

la figure 6 montre revolution (OF) de la fonction objectif en fonction du nombre k de 
simulation d'ecoulement executees, GDMl correspondant a une optimisation menee en 
combinant une realisation de depart et une seule realisation complementaire, 
GDMGBC3, a une optimisation menee en combinant la realisation de depart et une 
realisation composite construite & partir de trois realisations complementaires, 
GDMGBCIO, a I'optimisation menee en combinant la realisation de depart et une 
realisation composite construite a partir de dix realisations complementaires. et 
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GDMGBC30, a une optimisation meriee en combinant la realisation de depart et une 



Description detaillee de la methode 

La methode selon I'invention permet d'orienter a chaque iteration la construction 
de la chatne de realisations afin d'atteindi"e une direction de recherche souhaitable. La 
technique retenue tire avantage de I'information apportee par la direction de descente 
definie par les gradients de la fonction objectif /. Elle peut etre mise en oeuvre au moyen 
d'un simulateur numerique d'un type connu des gens de Tart tel que les simulateurs 
ATHOS ou ECLIPSE. 

Recherche orientee par la direction de descente 

On considere a ce stade la direction de descente (evolution vers une -^valeur 
minimale qu'elle soit ou non locale) definie par les gradients de la fonction objectif J par 
rapport a toutes les composantes de la realisation z- Ces gradients se deduisent des 
coefficients de sensibilite : 



Le probleme du calcul des coefficients de sensibilite a ete largement aborde dans la 
litterature scientifique. On pouira, par exemple, consulter : 

Sun, N.Z., Inverse problems in groundwater modelling, Kluwer Acad. PubL, 
Dordrecht, The Netherlands, 1994. - 

En particulier, la technique de J'etat adjoint permet de calculer Tensemble de ces 
coefficients en resolvant le probleme direct et son probleme adjoint, dans un temps 
equivalent a deux fois celui requis pour le resolution du probleme direct. 

Comme Tobjectif est d'integrer I'infonTiation foumie par les gradients dans la 
methode de deformation graduelle, il faut se ramener a la realisation Gaussienne qui resulte 
de la transformation de On distingue deux cas selon 1) que Ton dispose ou non de 
donn6es ponctuelles de conditionnement, c'est-a-dire de mesures de z en certains points. 



realisation composite construite a partir de trente realisations complementaires. 
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Dans le cas ou I'on ne dispose d'aucune donnee, si y est la realisation Gaussienne 
obtenue en transformant alors : 



dJ _ dJ 



Dans le cas contraire, on suppose que z est connu en certains points a. Ainsi la 
realisation Gaussienne s, deduite de la transformation de z, est une realisation 
conditionnelle de F. Pour obtenir s, on genere une realisation non conditionnelle y dc Yd 
on la conditionne aux valeurs connues en a par krigeage. La realisation conditionnelle se 
deduit de : 

S = 5*H-(y — y-^) 

5* et sont calculees respectivement par krigeage a partir des donnees reelles et 
les donnees de y generees au niveau des points a. On montre alors que : 



dJ 



0, Mn-^a 



Pour des proprietes physiques continues, I dy^^ ou dzy^ Ids^^ se calculent a partir 
de la fonction anamorphose qui peiTnet de transformer la realisation z en realisation 
Gaussienne. Lorsque les proprietes physiques considerees sont categorielles ou discretes, 
ces derivees n'existent pas. Les techniques de gradients ne peuvent alors s'appliquer. 

Ces differentes relations sont d'un interet direct si on se contente de pouvoir 
deformer globalement une realisation. Si on souhaite par contre la deformer localement, il 
faut appliquer la methode de deformation graduelle au bruit blanc Gaussien x utilise pour 
generer y, Dans ce cas, une etape supplementaire s'impose : le calcul des derivees de la 
fonction objectif par rapport aux composantes du bruit blanc Gaussien. 

Pour illustrer le calcul de ces derivees, on propose de se concenti^er sur le cas 
particulier suivant lequel la realisation. Gaussienne non conditionnelle y est obtenue a partir 
du generateur FFT-MA decrit dans Tarticle public par Le Ravalec et aL 2000 cite plus haut. 
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Le principe de base de ce generateur FFT-MA (pour "FFT-Moving Average") est 
de transformer un bruit blanc Gaussien x en realisation Gaussienne y correlee a partir d*un 
produit de convolution : 

5 La fonction g resulte de la decomposition de la fonction de covariance C telle que 

C = ^ * ^' , ou est la iransposee de g. Les derivees de la fonction objectif par rapport 
aux composantes du bruit blanc Gaussien sont : 

dJ dJ 3y, 

L' expression discrete du produit de convolution conduit a yi - ^^gi^k^k > Qi^i 

jt 

10 implique dyjdx^ =gi-n • Si on introduit cette formulation dans les derivees de la fonction 
objectif, on montre : 




Cette formule exprime le fait que la derivee de la fonction objectif par rapport au 
n*^™ composant du bruit blanc Gaussien est donnee par la /z^*^™ composante du champ 
15 obtenu en convoluant toutes les derivees de la fonction objectif par rapport aux 
composantes de la realisation Gaussienne avec le kernel de la fonction de covaiiance. A 
partir du cadre mis en place pour le generateur FFT-MA, on procede comme suit pour 
determiner ces derivees, 

1- calcul de la transformee de Fourier de dJ/dy, ces derivees etant foumies par 
20 mise en oeuvre du simulateur numerique direct ; 

2- multiplication de cette transformee de Fourier avec celle de g qui est foumit par 

FFT-MA ; 



3-calcul de la transformee de Fourier inverse du produit precedent. 
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Le temps requis pour le calcul de ces derivees est negligeable : i] represente un 
temps additionnel de 2/3 par rapport a la simulation d'une realisation Gaussienne par FFT- 
MA. 

Quels que soient le generateur de realisations et le simulateur numerique direct, on 
5 suppose, dans la suite, que Ton est capable de definir une direction de descente a paitir des 
derivees de la fonction objectif. Si Toptimisation de la fonction objectif est menee par 
rapport a cette direction de descente seulement, on detniit en general la coherence de la 
realisation vis a vis du modele de variabilite spatiale. Dans la section qui suit, on integre 
r information apportee par les derivees de la fonction objectif dans le schema de 
10 combinaison graduelle. 

Prise en compte des derivees de la fonction objectif dans le processus de 
deformation graduelle 

On considere la chaine de realisations y\{t) construite a paitir de yo d'une autre 
realisation w de 7 (figure la). A present, au lieu d'utiliser une realisation complementaire u 
15 telle qu'elle, on tire au hasard P realisations complementaires Up (p=l,2,.,.,P) de Y et on 
substitue a u une combinaison des Mp=i,p (figure Ic), Cette combinaison n'est pas 
quelconque : elle suit le mode de construction suivant : 

p p " 

La realisation u resultante est une realisation de Yei est independante de vo. D'apres 
20 Tequation 1, u est aussi la direction de recherche calculee pour la chaine yi(t) au point de 
depart yo. On propose done de construire u de sorte que la direction u soit aussi proche que 
possible de la direction de descente donnee par les gradients en yo. 

On definit au prealable I'espace V des vecteurs a N dimensions muni du produit 
scalaire ; 

25 {ynyj) = jTf,y,nyj.„ ^y^yj^v 

yi,n sont respectivement les n^^'''^^' composantes des vecteurs y^ et yj . Soit 

U un sous-espace de y defini par la base orthonormale (wj • La direction de 
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recherche. dans U la plus proche de la direction de descente v est donnee par la projection 
de vdans U (figure 3) : 

p 

En normalisant ce vecteur, on obtient la direction u souhaitee. Les coefficients de 
5 combinaison X de I'equation (2) sont done : 




La realisation u ainsi definie est appelee realisation connposite. 

On a considere jusqu'ici la construction de chaines de realisations par combinaison 
la realisation initiale avec une realisation composite. Cette technique peut toutefois etre 
10 generalisee a la constmction de chames faisant intervenir un certain nombre M de 
realisations composites, toutes obtenues par composition a partir de P;„ realisations de F, 

ces realisations etant independantes. L'optimisation fait alors intervenir M 

parametres. Cette technique augmente le degie de liberte dans le processus d'optirnisation, 

* - 

mais oblige a gerer M parametres d'optirnisation. 
15 Exemple numerique 

On constiiiit un modele synthetique de reservoir sur lequel on teste la'methode 
selon I'invention. 

Le reservoir synthetique de reference est represente sur la figure 4A. II s'agit d'un 
reservoir mono-couche comprenant 51x51 mailles de 10 m d'epaisseur et de 40 m de cote. 

20 La distribution de permeabilite est lognormaJe de moyenne 200 mD et d'ecart-type 100 
mD. Le logarithme du champ de permeabilite est caracterise ' par un variogramme 
spherique isotrope et une longueur de correlation de 480 m. Les autres proprietes 
petrophysiques sont constantes : la porosite est de 25 %, la compressibilite totale de S.IO"'^ 
bar'^ et la viscosite du fluide de 1 cP. Un puit de production BHP-PROl, de rayon 7.85 cm 

25 depourvu d'effet parietal ("sk.in effect'*) se trouve au centre du reservoir : il est entoure par 
quatre puits d'observation (BHP-OBSl, BHP-OBS2, BHP-OBS3, BHP-OBS4) (figures 4). Une 
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simulation numerique d'ecoulement permet d'obtenir pour ce resen^oir un jeu de donnees 
de reference comprenant les pressions aux cinq puits (figure 5). 

L'objet du probleme inverse est de deteraiiner un modele de reservoir coherent 
avec les donnees de pression en supposant inconnue la distribution des permeabilites. Dans 

5 ce but, on lance quatre processus d' optimisation en partant de la meme realisation de 
depart (figure 4B). Pour chacun de ces processus, on considere un unique parametre 
d'optimisation, a savoir le parametre de deformation. Le premier processus (GDMl) 
reprend le schema de deformation graduelle classique avec la construction d'une chaane de 
realisation en utilisant la realisation de depart et une seule realisation complementaire. Les 

10 trois autres processus (GDMGBC3, GDMGBCIO et GDMGBC30) illustrent rapplication 
de la methode selon Tinvention : les chaines sont dans ce cas elaborees avec la realisation 
de depart et une realisation composite issue de la combinaison de 3, 10 et 30 realisations 
complementaires. La realisation composite est construite comme explique ci-dessus (Eq. 2) 
en integrant I'information foumie par les gradients de la fonction objectif par rapport au 

15 bruit blanc Gaussien. Pour chacun des processus, on represente sur la figure 6 revolution 
de la fonction objectif (OF) en fonction du nombre k de simulation d'ecoulement 
executees. . . 



20 



On observe qu'en utilisant les gradients et en augmentant le nombre de realisations 
complementaires, la fonction objectif diminue plus yite pour un meme nombre de 
simulations d'ecoulements. 
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RE\^NDICATIONS 

1) Methode pour former plus rapidement un modele numerique stochastique de type 
Gaussien ou apparente, representatif de la distribution d'une grandeur physique dans un 

5 milieu heterogene poreux, qui soit cale par rapport a des donnees dynamiques obtenues par 
des mesures effectudes dans le milieu ou des observations prealables, et caracteristiques du 
. deplacement des fluides dans celui-ci, comportant un processus iteratif de deformation 
graduelle ou Ton combine lineairement a chaque iteration, une realisation initiale (y) d'au 
moins une partie du milieu, et au moins une deuxieme realisation independante de la 

10 realisation initiale, les coefficients (Oi) de cette combinaison lineaire etant tels que la 
somme de leurs canes est egale a 1; et Ton minimise une fonction objectif (J) mesurant 
recall entre un jeu de donnees non lineaires deduites de la dite combinaison au moyen d*un 
simulateur d'ecoulement dans le milieu, et les dites donnees geologiques et dynamiques 
par ajustement des coefficients de la combinaison, le processus iteratif etant repete jusqu'a 

15 obtenir une realisation optimale du modele stochastique, caracterisee en ce que Ton 
accelere la vitesse de deformation graduelle vers le modele optimal representatif du milieu, 
en selectionnant comme deuxieme realisation a combiner avec la realisation initiale, au 
moins une realisation composite obtenue en selectionnant au prealable une direction de 
descente definie en fonction des gradients de la fonction objectif par rapport a toutes^les 

20 composantes de la dite realisation initiale (y), 

2) Methode selon la revendication 1, caracterisee en ce que Ton obtient la 
realisation composite par combinaison lineaire d'un jeu de P realisations independantes du 
modele, les coefficients de la combinaison etant calcules de fafon a ce que la direction de 
descente a paitir de la realisation initiale y soit le plus proche possible de celle definie par 

25 les gradients de la fonction objectif par rapport a toutes les composantes de la realisation 
initiale. 

3) Methode selon la revendication 1 ou 2, caracterisee en ce que Ton conduit 
I'optimisation a paitir d'un parametre de deformation, qui controle la combinaison entre la 
realisation initiale et la realisation composite. 

30 4) Methode selon la revendication 1 ou 2, caracterisee en ce que la dite 

combinaison n'affectant qu'une partie de la realisation initiale, on applique le processus 
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iteratif de deformation graduelle a un bruit blanc Gaussien utilise pour generer une 
realisation gaussienne et on detennine les derivees de la fonction objectif par rapport aux 
composantes du bruit blanc Gaussien. 

5) Methode selon la revendication 1 ou 2, caracterisee en ce que Ton combine la 
realisation initiale avec un certain nombre M de realisations composites, toutes obtenues 
par composition a partir de Pm realisations independantes de Y, T optimisation faisant 
intervenir M parametres. 
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